Finger behavior of a shear thinning fluid in a Hele-Shaw cell 
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We make a theoretical study of the behavior of a simple fluid displacing a shear thinning fluid 
confined in a Hele-Shaw cell. To study the Saff man- Taylor instability when the displaced fluid 
is non Newtonian we face the problem of having a field which is non laplacian. By means of an 
hodographic transformation we are able to solve the problem in the case of weak shear thinning 
while taking into account the non laplacian character of the equation. Our results predict that the 
finger width decreases towards zero for small values of the surface tension parameter which inversely 
proportional to the finger velocity. 

PACS numbers: 47.20.-k, 68.f0.-m, 83.10.Lk 
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Instabilities in complex fluids have been recently the 
subject of intense research due to their importance in 
technology. The classical hydrodynamic experiments, 
such as the Dean flow and Taylor-Couette flow have 
been considered to study bulk instabilities On the 
other hand, interfacial instabilities between complex flu- 
ids, with some exceptions [0, have been less studied. 
Theoretically, this is due to the difficulty of the equations 
governing these systems. Here we consider the Saffman- 
Taylor problem for a viscous fluid whose viscosity changes 
with shear. The Saffman- Taylor problem H|J is the pre- 
diction of the steady-state shape of the fluid interface in a 
two-phase flow confined in a linear Hele-Shaw cell. A vis- 
cous fluid is pushed by a low viscosity fluid; modes grow 
and compete dynamically, and the competition leads to 
a single finger-shaped pattern at large times Q . Experi- 
mentally |^j6|, the finger is characterized by its width A. 
It is found that A > 1/2 in units of the channel width 
W, and that it is a unique function of a control param- 



eter a such that A 



as a — > 0. The surface-tension 



parameter a is defined as 



4T 
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where T is the surface tension, /j, is the viscosity of the 
fluid that is being pushed, b is the gap spacing of the 
Hele-Shaw cell and U is the velocity of the finger in the 
laboratory frame. The control parameter a provides a 
singular perturbation to the zero-surface-tension equa- 
tions of motion describing the Hele-Shaw flow Qj. Several 
experiments were carried out to alter this parameter and 
study the system's response |p|-^2[. New steady-state 
patterns were observed both in experiments and simu- 
lations in these "perturbed" Hele-Shaw cells 
Most of this patterns were theoretically explained, by 
altering the boundary conditions of the problem |j 14 
while keeping the Laplacian character of the equations 



(see below). Other patterns were explained by impos- 
ing mathematical conditions at the finger tip |fi5[ . Re- 
cently, experiments have been made in which the fluid 
which is pushed is non Newtonian. Some of these ex- 
periments have used clays [ fl6| and others polymer solu- 
tions 0. Both groups report seeing a fracture- like be- 
havior of the complex fluid. Non-Newtonian fluids differ 
widely in their physical properties, with different fluids 
exhibiting a range of different effects from plasticity and 
elasticity to shear thickening and shear thinning. In this 
paper we are interested in studying how the viscous fin- 
gering regime is altered when considering an effect such 
as shear thinning or shear thickening, i.e., we consider 
a non-Newtonian fluid whose viscosity depends on the 
local shearing. 

Classically the constant- velocity displacement of a vis- 
cous fluid by a fluid of negligible viscosity such as air is 
governed by Darcy's law, which states that the velocity 
is proportional to the pressure gradient. Recently a gen- 
eralized Darcy's law where viscosity depends upon the 
squared pressure gradient has been derived by using a 
fluid model with shear rate dependent viscosity jL8) . As 
pointed out by them, this can be formally inverted to 
give a generalized Darcy's law where viscosity depends 
on the magnitude of the local velocity. We consider the 
latter for the present work. For a constant viscosity fluid, 
Darcy's law plus the incompressibility of fluids lead to a 
Laplace's equation for the pressure. This allows one to 
use well established conformal mapping techniques [fl9|| . 
For a non-constant viscosity fluid, imposing the incom- 
pressibility of fluids does not lead to Laplace's equation 
for the pressure. This poses a serious problem since one is 
unable to apply standard techniques. In the laboratory 
frame, where v v i a t es = 0, the equations governing our 
system are the generalized Darcy's law where viscosity 
depends on the magnitude of the local two-dimensional 
fluid velocity and the incompressibility of fluids: 
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With the following boundary conditions 



V -v = 



v n = UsinOr , r 

T > on tne linger 



(2) 



(3) 



and 



v y = on the walls, v. 



— UX, v y — as x — > oo 
for —1 < y < 1 Uj. = 0, = a; ^ — oo for A < 
\y\ < 1. Here i? is the local radius of curvature, V is the 
magnitude of the local fluid velocity in the laboratory 
frame of reference, Or is the angle tangent to the finger 
measured from the x-axis and the subfix n stands for the 
component of the velocity normal to the interface. We 
have set W = 2. It is worth it to explicitly point out that 
equation |^ does not violate Galilean invariance since it 
is only valid in the laboratory frame of reference. We 
are going to treat problems for which viscosity depends 
on velocity, so we write the viscosity as ^ = -j^(j^), 
where p depends on the velocity and po is the function 
p evaluated at the finger tip, po is a constant defined in 
such a way as to be the viscosity of the fluid when there 

P, 
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is no dependence on the velocity. Defining 
equations || can be written in a differential form as 

Vydy) , (4) 



— (v x dx 
P 



and 



dx 



-Vydx 



v x dy 



(5) 



We write the complex velocity in polar coordinates as 
(v x — iv y ) = V exp —id and express dz = dx + idy as 



dz = — ^— (— deb + idx) 
V po 



(6) 



The basic problem that we face is the non-Laplacian char- 
acter of the equation describing the problem, i.e., equa- 
tions ||do not lead to Laplace's equation for the pressure. 

We make an hodographic transformation by exchang- 
ing dependent and independent variables. So we consider 
x(V, 9) and y(V, 6). We express d<b and d\ in terms of 
partial derivatives as functions of (V,9). Since dz is an 
exact differential its crossed derivatives should be equal. 
This gives an equation whose real and imaginary parts 
provide identities among the partial derivatives of x an d 
<b with respect to V and 6. Such identities allow us to 
write equations for x(Vi @) an( l <t>(Y, 0) by equating the 
crossed derivatives obtained by differentiating them. 

We work out the case in which the viscosity depends 
on the velocity as a power law, i.e., p — V a . This is one 
of the standard forms for the viscosity found in literature. 
In this case, the equations for xiYi 6) an d ^(V, 6) have 
the following form 



d 2 x (1 - a) dx 1 d 2 x 



dV 2 



V dV V 2 dr 2 
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dV 2 



V dV 



±^± = 
V 2 dr 2 



(8) 



Where the angle has been redefined as r = (1 — a) -1 / 2 ^. 
We rescale the fields x an d (b in such a way as to obtain 
Laplace equation to linear order in a, i.e., for the fields 



V = V a/2 (b and X = V^" /2 
we have the following equations 



X 



d 2 X 1 dX 1 d 2 X a 2 X 



and 



dV 2 V dV V 2 dr 2 4 V 2 



d 2 V 1 dV 1 d 2 V _ a 2 V 
W 2 + VdV + V 2 ^ 2 ~ IF 
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Up to this point we have shown that, for the viscosity 
depending on a power of the velocity 



V 2 V V = 0(a 2 ) and V 2 V X = 0(a 2 ) 



(12) 



where Vy is the Laplacian operator in the velocity plane. 

The hodograph method and the rescaling of fields, de- 
spite being little intuitive from a physical point of view, 
have permitted us to write the equations of the problem 
in such a way as to allow for an approximation in a con- 
trolled manner. Now, for small a which corresponds to 
the case of weak shear thinning, we approximate V and X 
by laplacian fields. This leads to the following expression 
for dz 



exp zt(1 — %) , 

dz = VU) ^ 



idX) 



(13) 



When writing the differentials for V and X from [13| it 
is easy to see that, if we define an effective velocity of 
the form v* = v*i + v*j with v* — y( 1- t) cosr(l - ^) 
and v* — I/( 1_ f ) sinr(l — |r) we have equations that are 
mathematically equivalent to classical Darcy's law and 
the incompressibility of fluids, i.e., 
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We emphasize here that v* is not a real velocity; and that 
V is not a pressure, but a function of pressure and the 
local fluid velocity. Nothing new has been introduced by 
defining v* . We have found equations which are mathe- 
matically equivalent to the classical problem, and there- 
fore, we are now dealing with Laplacian fields in real 
space. We are now in the position of performing confor- 
mal transformations. We work in a frame of reference in 
which 



^ v-^n , T x + mi 

$ = and * = 



(15) 
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where 7i is an unknown analytic function whose val- 
ues at the contour are set so that the contour 
(wall+finger+center line) maps into the real axis of the 
(s,t) plane of the following conformal transformation 



= s + it = cxp - $ + i* + iT~t]n 



(16) 



This method has been used previously in reference P0| . 
McLean and Saffman have solved the problem of constant 
viscosity for arbitrary surface tension . In the rest of 
the analytical treatment we follow closely their technique. 
Details of the derivation will be presented in a longer 
publication. In the variables 9 = 9t—tt and = <z(l — A) 
we find 



ksqn 



Qn 



ds 
q N U 
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where 



1-A ttQo/AJo Qn(s'~s) 



Tb 2 ir 2 



qN 



-ds' 
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q = ^ where dS means differentiation along the inter- 
face and Qo = V^ a ^ 2 . Equation [l7| reduces to McLean- 
Saffman equation when a = 0. A second equation is 
obtained by analyticity properties and it is exactly like 
the McLean and Saffman relation among q^ and 9, i.e., 



log q N (s) 



6(s') 



-ds 1 



(19) 



io s'(s'-s) 

Note that in this case qN no longer represents the finger 
velocity at the finger-tip frame. The third equation is 
obtained by writing explicitly the fact that the square 
of the velocity is equal to the sum of the squares of the 
normal and tangential components, and it is 



KsqN 



yl-a/2 



ds 
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sq N 
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FIG. 1. Finger width as a function of the surface tension 
parameter k. The solid line is the line of constant viscos- 
ity (a = 0), the other lines correspond to the case of shear 
thinning with exponents: o a = 0.2, + a = 0.4, n a = 0.6. 

Equations [Tt], [jl], and generalize the Saffman- Taylor 
problem to the case in which viscosity depends as a power 
law on the local velocity. They are solved numerically 
with the method described in fll9[| . In figure [l] we plot 
finger width as a function of k for the case of shear thin- 
ning (a > 0) which is the case more often seen in ex- 
periments. We also plot for reference the case a = 0, for 
which we obviously recover the McLean-Saffman solution 
to the problem since the viscosity is a constant. At large 
values of k the solutions are very close to the solutions 
of McLean-Saffman. On the other hand, at small values 
of k the solutions decrease dramatically towards a zero 
finger width. Between these two regions there is an in- 
termediate region in which solutions go smoothly from 
one type of behaviour to the other one. In all cases, solu- 
tions deviate from the a = solution and the deviation 
is stronger the larger the value of a. We present results 
for a = 0.2,0.4,0.6 in order to illustrate how the solu- 
tion changes with the power a. We have also performed 
a solvability analysis in the case of shear thinning. Such 
analysis gives a scaling among the finger width A and the 
surface tension a in the small a limit. The result is 



A 



*/((!+«) (2-a)) 



(21) 



It is unfortunate that numerically it is impossible to ver- 
ify such a scaling, since it would imply computations at 
values of a out of our numerical capabilities. Therefore, 
this result should be taken as complementary to our nu- 
merics. 

We have treated the problem of viscous fingers in the 
presence of shear thinning. We have obtained fingers 
whose width decreases towards zero as the surface tension 
parameter goes to zero. A detailed derivation of the equa- 
tions presented here, comparison with the ?7-model |2l]] 
and a description of the solvability analysis will be pre- 
sented in a longer publication. In this forthcoming publi- 
cation we will present as well results for shear thickening 
which, surprisingly, present a dramatic decrease of finger 
width at very small n. Also, since our model is two di- 
mensional, it is valid as long as the finger width is larger 
than the height of the cell. Effects due to a third dimen- 
sion, which might be very important, specially in con- 
nection to fracture behavior, have not been considered in 
the present work. 
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